Quantum dynamics in two and three-dimensional quasiperiodic tilings 
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We investigate the properties of electronic states in two and three-dimensional quasiperiodic struc- 
tures : the generalized Rauzy tilings. Exact diagonalizations, limited to clusters with a few thou- 
sands sites, suggest that eigenstates are critical and more extended at the band edges than at the 
band center. These trends are clearly confirmed when we compute the spreading of energy-filtered 
wavepackets, using a new algorithm which allows to treat systems of about one million sites. The 
present approach to quantum dynamics, which gives also access to the low frequency conductivity, 
opens new perspectives in the analyzis of two and three-dimensional models. 
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Since the discovery of quasicrystals 0], the nature of 
the eigenstates and of the energy spectrum, as well as 
the propagation of electrons in quasiperiodic potentials, 
have been the subject of many theoretical studies. One- 
dimensional (1-D) quasiperiodic Hamiltonians have been 
widely investigated by means of numerical and per- 
turbative approaches QJ^. Their spectrum consists in a 
set of zero width bands associated to critical eigenstates 
responsible of a sub-ballistic propagation of the electrons. 
In higher dimensions, the situation is more complex and, 
given the geometrical complexity of the structures, ana- 
lytical treatments are difficult. Most of the studies have 
focussed on topologically trivial systems whose proper- 
ties can be easily extracted from the ID case (Fibonacci 
quasilattices j^J^, labyrinth tihng [^|j9|). It has been 
shown that if the propagation is always sub-ballistic, the 
spectrum could be either absolutely continuous, singular 
continuous or any mixture. Other works based on exact 
diagonalizations, have also been performed on Penrose- 
like 2D and 31? systems. In 2D, there is some evidence 
that eigenstates are critical [|o| and responsible of an 
anomalous diffusion [0. In 3D systems, the nature of 
eigenstates remains unclear fl^ . In this context, it is 
of great interest to determine whether critical states are 
generic of true quasiperiodic systems especially in the 
3-D case, and to study the physical properties induced by 
such states. 

In this letter we study in a tight-binding approach, 
the electronic properties of the of the 2D and 3-D gene- 
ralized Rauzy tilings (GRT) [Q. We discuss the nature 
of the eigenstates by computing their localization degree 
for small systems (~ 5500 sites). This analyzis suggests 
that the eigenstates are critical and rather more extended 
at the band edges than at the band center. To go beyond 
these exact diagonalizations, we introduce a powerful al- 
gorithm that allows us to study the quantum dynamics 
for very large systems (up to 10^ sites here). Sub-ballistic 
propagation laws are clearly established and the energy 



dependence of the exponents confirms that states are 
more mobile (almost ballistic) close to the band edges 
than in the center. From this point of view, these results 
provides evidence of critical states in the 2D and 3-D 
GRT. Finally, we discuss the low frequency conductivity 
which can be directly related to the quantum diffusion 

The GRT are codimension one quasiperiodic structures 
that can easily be built in any dimension D by the stan- 
dard cut and project method [|l3|. Here, we focus on 
approximant structures that are obtained for rational 
cut directions. These tilings have a complex topologi- 
cal structure with sites of coordination number ranging 
from D + 1 to 2D + 1. as displayed in Fig. ^ for D = 2. 




FIG. 1. A unit cell of the 10* 
with sites. 



approximant of the 2D GRT 



This geometrical complexity implies that they cannot 
be separated into a combination of 1-D systems as in the 
Fibonacci quasilattices or as in the labyrinth tiling 
Q. We thus expect that their spectral and dynamical 
properties are representative of other 2D or 3-D Penrose- 
like quasiperiodic structures. 
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We consider noninteracting electrons described by the 
following tight-binding Hamiltonian : 

where t = 1 is the hopping energy between nearest neigh- 
bor sites, and where \i) denotes an orbital localized on 
the site i. As discussed in the Hamiltonian matrix 
with periodic boundary conditions is straightforwardly 
written as a band (Toeplitz-like) matrix. This is due 
to a natural indexing of the sites known as the conum- 
bering which classify them according to their local envi- 
ronment. The study of the spectrum of H has already 
been presented in Ref. | [T7[ | , and we focus here on the na- 
ture of the eigenstates. Therefore, we have computed by 
exact numerical diagonalization, the participation num- 
ber P of each normalized eigenstate Itp) defined by : 

= (X^i I I'*) • For a given system size, this 
quantity measures the localization degree of the state {ip) 
considered. We display in Fig. (||), for the 2D and 'SD 
GRT, the participation ratio p = P/N as a function of 
the energy (N is the number of sites per unit cell) . Qua- 
litatively, we observe in both systems that, despite strong 
fluctuations, the most extended eigenstates are globally 
located near the spectrum edges and the more localized 
ones in the central part. 
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FIG. 2. Participation ratio as a function of the energy for 
the 14*'' approximant (5768 sites) of the 2D GRT (left) and 
for the 13"' approximant (5536 sites) of the 3D GRT (right). 

Another interesting information is given by the sca- 
ling of the mean participation number (averaged over the 
whole spectrum) versus the number of sites N. Indeed, 
when the system size increases, the mean participation 
number converges towards a finite value if the states are 
exponentially localized and goes to infinity otherwise. As 
shown in Fig. ^, this mean participation number behaves 
as N'f with ~ 0.965 ±0.001 and 73D = 0.975 ±0.003, 
indicating a weakly critical behaviour of the eigenstates. 
Note that this is only an average behaviour that conceals 



strong fluctuations from one energy to another as it has 
already been shown in Penrose- like tilings . 

In the 2D GRT, the exponent is greater than the 
one obtained for the octagonal tiling ||l8| (codimension 
2) where 7 = 0.87 ± 0.05. This difference can be un- 
derstood by invoking the codimension of the structures 
that gives an indication on the geometrical complexity 
in terms of local environment. Intuitively, we indeed ex- 
pect low codimension structures to be weakly quasiperi- 
odic, and to have weakly critical states. The recent level 
statistics analyzis done in the 2D GRT |jl^ also corro- 
borates this assumption. Finally, up to the numerical 
precision, we have "f2D < "fsD- It means that for a fixed 
codimension (here 1), the localization of the eigenstates 
increases when the dimension is lowered, as it is the case 
in disordered systems. To confirm the trends observed 
with diagonalizations, we have studied the spreading of 
energy-filtered states on much larger systems of about 
one million sites. 
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FIG. 3. Mean participation number as a function of the 
number of sites for the 2D (□) and the 3D (O) GRT. 

The diffusion of wave packets can be measured by the 
mean square spreading : 

X^E,t)^{iXit)-X{0)f)E (2) 
= (AHt)Ait))E, (3) 

where A{t) = [X,e~^^'^]. X{t) is the Heisenberg repre- 
sentation of the position operator along the x direction 
and (...)£; denotes the average value over the eigenstates 
of energy E. We stress that X'^{E,t) evolves with time 
even though it is a mean value over eigenstates. 

X^{E, t) is related to (t{Ef, oj) the conductivity at fre- 
quency Lo and Fermi energy E-p [ p^ . In the limit of small 
frequencies one has : 

{a{E,,u;))^^e'l e^^' (^niE^ f ^dt^ ''^ ) ^4) 
^ 

where n^Ep) is the density of states at the Fermi en- 
ergy, e is the electron charge, and (. . .)a means an av- 
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erage over a range A of values of Ep. Note that at zero 
frequency, the integral in is equal to the diffusiv- 
ity D{Ep) = liiiit^ocdX^{EF,t)/dt and one recovers the 
Einstein formula a{Ep,0) = e'^n{Ep)D{Ep). 

To investigate the quantum diffusion, we consider an 
approximant of the GRT with a sufficiently large unit 
cell so that the eigenstates are Bloch waves characterized 
by a wave vector k. As for the exact diagonalizations 
we restrict our study to the k = subspace, i.e. to the 
eigenstates which are invariant by translation from one 
unit cell to another. Note that for an approximant the 
operator A{t) is invariant by translation of a unit cell 
and thus lets invariant the subspace of k = states. We 
calculate X'^{E,t) defined by : 



Xl{E,t) = {A\t)A{t))EM=o 

Tv{A'^{t)5{E - H)A{t)] 



X^iE,t) 



Tt{S{E- H)} 



(5) 
(6) 



Actually, Xq{E, t) is a good approximation of X^{E, t) 
provided that t is lower than the time required to signif- 
icantly couple orbitals separated by a distance lx/2 dur- 
ing the time evolution. In this study, we have always 
considered wave packet evolutions for such times. The 
operator trace Tr is restricted to the subspace k = 0. 
An efficient trick to evaluate the traces is to use states 
\ip) having a constant presence probability on each site 
I (il^/i) |2 = 1 /jv and a random phase on each orbital of the 
unit cell. We thus approximate Ti{B} by N{ip\B\^). In 
fact, we should perform an average over several random 
phases configurations but we have checked that, thanks 
to the great number of sites in the unit cell, we only need 
to consider one random phase state \ip) to get a good es- 
timate of the traces. Moreover, to eliminate fluctuations, 
we have performed a convolution of the 6{E — H) by a 
Gaussian operator of typical width A = VF/100 where W 
is the total bandwidth This is equivalent to com- 

pute the spreading of energy filtered wavepackets. These 
parts of the calculation are rather standard and in fact 
the crucial step is to evaluate A{t)\^p). 

We start from the expansion of the evolution operator 
Q-iHt ygjj^g ,^ basis of Chebyshev polynomial Pi{E) of 
the first kind. These polynomials are associated to the 

This expansion 



weight p{E) = 1/ (^7r^462 - {E ~ af ' 

converges provided the spectrum of H (computed here 
using the standard recursion method [pT[) is included in 
the interval [a — 2b,a + 26] ||2^ . The decomposition of 
the evolution operator is written as : 



-IHt 



oc 

E 

i=0 



Cr{t)P,{H), 



(7) 



with Cp{t) = i^hpJp{—2bt) where ho — 1, hp^o = 1 /2 , and 
where Jp{x) is the Bessel function of order p |^,^. At 
large p and fixed x, one has : Jp{x) ~ l/y^2np{xe/2p)P 



so that the sum in Eq. (0) converges quickly and can 
be truncated. The polynomials Pi obey the following 
relation : 



bP,+i{H) = {H- a)PdH) - bP,^i{H) 



(8) 



with Po{H) = 1 and Pi{H) ^ {H - a) /2b. Then, intro- 
ducing 10,) = P^{H)\i;) and It/-,) = [X,P,{H)]\i^) and 
using Eqs. (Q) and (||), one readily obtains : 

b\4>,+i) = {H-a)\dp,)-b\^,-i) (9) 
h\ij,+i) = {H- a)\^i) - felV',-!) + [X,H]\^i). (10) 

Solving this latter system allows one to compute efh- 
ciently : 



A(OI^) =^C.(t)|Vz), 



(11) 



i=0 



and thus (E., t) with an arbitrary accuracy given by 
the order of the truncation. 

All the results discussed here have been obtained in the 
22"'' approximant (755476 sites) for the 2D GRT and in 
the 21"* approximant (1055026 sites) for the iD GRT. We 
have computed, for both tilings and for several energy- 
filtered wave packets, the time dependence of Rq{E, t) — 
X^{E, t) -t- Y§{E, t) + Z^{E, t) where a;, y and z are three 
orthogonal directions {Zq = in 2D). 
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time time 
FIG. 4. Behaviour of R^{E,t) for three energy-filtered 
wave packets in the 2D GRT (left) and in the 3D GRT (right). 

As displayed in Fig. |, one has Rl{E,t) ~ t2/3(£;) 
large times with < /?(£') < 1. The observation of 
such a sub-ballistic propagation for large approximant is 
consistent with the existence of critical eigenstates. The 
energy dependence of the diffusion exponent (3 shown in 
Fig. ^ also clearly confirms the tendency observed with 
exact diagonalizations since the diffusion exponent is, in 
both systems {2D and 3D), larger at the band edges 
than at the spectrum center. In the iD GRT, we have 
even obtained an almost ballistic propagation {(3 ~ 1). 
In addition, the shape of the density of states is nearly 
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free- electron like at band edges which is also consistent 
with the less localized eigenstates at these energies. If 
this general trend is similar to what is observed in the 
octagonal tiling and in the Penrose tiling , it is the 
opposite of what occurs in the icosahedral Amman- 
Kramer tiling [Q. Note also that in disordered systems, 
states at band edges are more localized than in the band 
center. The parameters that determine the nature of 
band edges states are obviously not yet understood. 



systems. Note that in the context of large scale elec- 
tronic structure calculations, order N algorithms raise 
much interest and can also be much more efficient than 
exact diagonalizations |Q. As a consequence, this ap- 
proach should open new perspectives in the investigation 
of 2D and 3D models. 

We thank Claude Aslangul for a careful reading of the 
manuscript. 
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FIG. 5. Normalized density of states (lower curve) and dif- 
fusion exponents (upper curve) of the 2D GRT (left) and of 
the 3D GRT (right). 

We would like to emphasize that the scaling laws are 
vahd for Rq{E, t) varying by 3 — 4 orders of magnitude in 
2D and 2—3 orders of magnitude in 3D. In both cases the 
number of sites N contained in a volume Rl{E,t) in 2D 
and Rq{E, t) in 3D varies by 3 — 4 orders of magnitude. 
This is much better than with exact diagonalization for 
which the typical variation of A'^ is 1 order of magnitude. 

The power-law behaviour obtained for RQ{E,t) allow 
us to explicitely calculate the low frequency conductiv- 
ity. Indeed, assuming X'^{EF,t) ~ C(£'F)i2/3(^^F), one 
obtains, using Eq. (^ : 

ct{Ef,lu) ~ e^niEp) C{Ef) T (2/3 -f 1) ii/uj)^^-\ (12) 

where F is the Euler Gamma function. For the GRT, 
the diffusion exponent P{Ey) is always greater than 1/2 
(superdiffusive propagation) so that although the system 
is non periodic, its dc conductivity a{EF,uj = 0) is infi- 
nite for all Ep. Note that experiments are qualitatively 
explained if /3 < 1/2 [|l6| but such values of /3 probably 
require a strong quasiperiodic potential. 

To conclude, we would like to point that the calcula- 
tion of quantum diffusion presented here, implies a com- 
putational load proportionnal to the number A^ of sites 
of the system whereas exact diagonalizations implies a 
load which is of order A^ or A^^ In A^ for the best algo- 
rithms. This is the main reason of the great efficiency 
of the present method and of its ability to treat large 
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